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Quantum tunneling often allows pathways to relaxation past energy barriers which are other- 
wise hard to overcome classically at low temperatures. However, this is not always the case. In 
this paper we provide exactly solvable examples where the barriers each system encounters on 
its approach to lower and lower energy states become increasingly wide and eventually scale 
with the system size. If the environment couples locally to the physical degrees of freedom 
in the system, tunnelling under these barriers requires processes whose order in perturbation 
theory is proportional to the width of the barrier. This results in quantum relaxation rates 
that are exponentially suppressed in system size; For these quantum systems, no physical 
bath can provide a mechanism for relaxation that is not dynamically arrested at low tem- 
peratures. The examples discussed here are drawn from three dimensional generalizations 
of Kitaev's toric code, originally devised in the context of topological quantum computing. 
They are devoid of any local order parameters or symmetry breaking and are examples of 
topological quantum glasses. We construct systems that have slow dynamics similar to either 
strong or fragile glasses. The example with fragile-like relaxation is interesting in that the 
topological defects are neither open strings or regular open membranes, but fractal objects 
with dimension d* = In 3/ In 2. 

Keywords: glassiness; stochastic processes; quantum tunneling 

1. Introduction: How slow is slow? 

The problem of the approach to equilibrium is a formidable one. Nature provides 
us with plenty of examples of systems that simply do not equilibrate in any exper- 
imentally accessible times, starting from a material that is so commonly accessible 
to become synonym of an entire field of research: Si02, glass [1]. A great deal of 
effort has been directed, both by physicists and chemists alike, to the question of 
how we can understand and describe these systems when they encounter dynamical 
obstructions in their attempted path to thermal equilibration. In spite of much re- 
search, however, a complete theoretical description of the glass transition remains 
an open problem. 

Understanding the origin of long equilibration time scales is both of fundamen- 
tal significance from a theoretical point of view and of great importance from a 
technological point of view. Uncovering the reasons why systems fall out of equilib- 
rium might enable us to manipulate them, either hindering or enhancing the slow 
dynamics depending on the specific application. Advancements on understandy- 
ing the properties of glassy materials would have implications to problems ranging 
from biology to petroleum recovery. 
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Before we begin our discussion, one should decide on a working definition of 
slow dynamics: how slow is slow? On the one hand, one can set a threshold - 
for instance, the famed viscocity of 10^^ Poise [l| - beyond which a system is 
said to be in a glassy state. A classification can then be developed based on how 
quickly the time scales grow as a function of some natural tuning parameter that 
leads the system into such state (e.g., temperature or concentration). According 
to this classification, if the parameter is temperature, one distinguishes between 
strong and fragile glasses if the equilibration time scales grow in an Arrhenius form 
T ~ exp[A/T] or faster, respectively. Examples of the latter are the exponential 
inverse temperature square behavior, r ~ exp[A^/T^], and the Vogel-Pulcher law 
for a true glass transition at finite Tq, t ~ exp[A/(T — Tq)]. 

On the other hand, one can classify systems in their glassy state by looking at how 
fast their large time scales grow with system size. While in a finite system all time 
scales are strictly speaking finite0 they can exhibit a dependence on system size. 
In general, one can account for three different scenarios, (i) There could be no size 
dependence at all. For instance, a purely Arrhenius behaviour with a local energy 
barrier A can lead to time scales ~ exp(A/r) that are independent of the size of 
the system, (ii) Time scales can grow polynomially in system size, r ~ poly(L). For 
example, this is the case of critical slowing down and diffusive modes. Finally, (iii) 
time scales can grow exponentially in system size, r ~ exp(L)H [In the following 
we shall use the short hand notation r ~ exp(L) to denote generic exponential 
time scales of the form exp(aL^.] Examples include the Sherrington-Kirkpatrick 
model [1] and p-spin glasses [1, 0] . 

Drawing a line between time scales that grow polynomially vs. exponentially in 
system size can be understood within the framework of computational complexity: 
Nature itself is a computer (which could in turn be simulated in a universal digital 
computer), running its algorithm for dynamical evolution. Any system thermalizes 
if the algorithm is allowed to run indefinitely. However, in finite times one has 
to come to terms with how efficient the dynamical evolution algorithm is. The 
system size L is the size of the problem and systems whose equilibration times 
grow exponentially in L correspond to computationally hard problems for nature's 
algorithm. 

In this article we focus on the system size dependence of equilibration time 
scales, in particular in quantum systems. In Sec. O we argue that any classical 
system whose relaxation is exponential in L, when endowed with local quantum 
dynamics, enters a quantum glass state (exponential in system size). The converse 
is not true. Quantum glasses exist that 'melt' at any finite temperature, with the 
appearance of time scales that are polynomial in system size. In Sec. [3] we provide 
three examples of this behaviour. 

We note that exponentially long times naturally arise in systems undergoing 
spontaneous symmetry breaking. For example, going from positive to negative mag- 
netization in the ordered phase of an Ising ferromagnet takes an exponentially large 
time in system size. These states however are distinguished by local order parame- 
ters (e.g., the local magnetization), which allow one to make significant progress in 
understanding their dynamics in terms of nucleation, domain growth and coarsen- 
ing. In this paper we focus instead on systems where exponential relaxation time 
scales appear in the absence of local order parameters and symmetry breaking. At 
the quantum mechanical level, a natural context where to look for ground states 



^ Note however that these time scales may already be in practice longer than any experimentally accessible 
times. 

^In any realistic (i.e., extensive) system, time scales cannot grow faster than exponential in the volume of 
the system. 
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without local order and without explicit disorder in the Hamiltonian or bath pa- 
rameters is that of topologically ordered systems jB, 0|. Indeed the examples of 
quantum glasses we put forward in Sec. [3] are inspired by toy lattice models for 
topological order [3] and ought perhaps to be called topological quantum glasses. A 
brief account on the physics of this family of quantum glasses, and inparticular of 
the example in Sec. 13. 2^ was discussed by one of the authors in Ref. [8(. 

Finally, in Sec. H] we briefly comment on general approaches to characterise dy- 
namical quantum systems and how - due to the characteristic exponential depen- 
dence on system size - one should beware of truncating perturbative approaches 
when studying quantum glasses, as they lead to reducible dynamics within discon- 
nected sectors. 



2. Classical and quantum glasses 

It is often the case that quantum tunneling allows relaxation past energy barriers 
which are otherwise hard to overcome classically at low temperatures. For instance, 
the growth of equilibration time scales due to Arrhenius activated behaviour from 
a local, finite barrier r ~ exp(A/r) is eventually cut off at sufficiently low temper- 
atures by temperature-independent tunneling across the barrier, Tg ~ A/t where t 
is some quantum mechanical amplitude for the process 

In general, classical systems that enter a glass state exhibit equilibration time 
scales that grow exponentially in system sizeH The origin of such characteristic 
exponential dependence lies in the appearance of large energy barriers between the 
glassy free energy minima, whose height and width grow with the system size. 

What happens if we take a classical system with r ~ exp(L) and lower the 
temperature to zero, where coherent quantum mechanical processes become active? 
As we shall argue hereafter, the fact that the barriers grow with system size has 
a dramatic effect on the quantum mechanical relaxation time scales. Realistically, 
we assume that only local terms are allowed in the tunneling Hamiltonian - for 
instance, a transverse field of magnitude t that flips individual spins in a localised 
spin system. Wide barriers mean that the system has to visit a large number 
Ng ~ poly(L) of excited states in its journey from one side to the other of the 
barrier. Quantum mechanically this process is strongly suppressed: Tq ~ {U/t)^% 
where U is the energy scale of the intermediate virtual states. Consequently, the 
quantum mechanical relaxation time scale acquires also a system size dependence 
of exponential form, Tq ~ exp[A'^s ln(C//t)] ~ exp(L). 

We conclude that a classical model with r ~ exp(L) remains exponentially slow 
at T = even if (local) quantum tunneling processes are allowed. A whole family 
of quantum glasses can thus be derived directly from classical glassy systems by 
replacing temperature with local quantum dynamics. One might wonder whether 
the converse is true: Are there quantum glasses that do not have a classical parent? 
In other words, are there quantum systems that have exponential relaxation times 
Tq ~ exp(L) at T = 0, but that immediately 'melt' to r < poly(L) as soon as 
T > 0? In Sec. [3] we will answer positively to this question by explicit construction 
of local Hamiltonians that are examples of purely quantum glasses. 



^ These time scales are dimensionless for convenience. They are intended as factors multiplying the charac- 
teristic microscopic time of the system in the absence of barriers. For instance, in the quantum mechanical 
tunnelling case, the characteristic time would be dictated by the inverse hopping amplitude, 1/t. 
^Systems where the glass transition occurs in the zero temperature limit ought to be treated with care, as 
the order of limits (thermodynamic vs. T — >■ 0) matters. For instance, time scales that are exponential in 
system size with a trivial Arrhenius activated prefactor exp(A/T) diverge in the T — >■ limit irrespective 
of system size. 
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3. Examples of quantum glasses without disorder 

In order to construct our quantum Hamiltonians, we shall draw inspiration from 
two rather distinct areas: (i) classical kinetically constrained models, where slow 
dynamics appear without disorder; and (ii) spin models for topological order, which 
exhibit gapped quantum ground states devoid of any local order parameters or 
symmetry breaking. 

Kinetically constrained models allowed to make progress on the study of the clas- 
sical glass transition by investigating systems in which the equilibrium properties 
are easy to understand, but the dynamics are non-trivial due to kinetic constraints 
(for a review, see Ref. [o^). Simple lattice models are constructed embodying the 
notion that there are jammed and unjammed regions in glass formers described 
by a discrete state, which is representable in terms of Ising spin variables [lo| . 
One particularly interesting class is that of spin plaquette models with nontrivial 
classical energy function and unconstrained spin dynamics, which can be equiva- 
lently described in terms of non-interacting defect variables with rather constrained 



multi-defect dynamics |lll4l4i|. The thermodynamics of these systems is therefore 
trivial, while their dynamics is not, displaying a rich non-equilibrium behaviour 
(ageing, for instance). 

Topological order is a relatively recent concept in systems of strongly interacting 
particles [1, [^. Some quantum phases of matter, in contrast to common examples 
like crystals and magnets, are not characterized by a local order parameter and 
broken symmetries. Instead, they are characterized by a ground state degeneracy 
(when the system is defined on a torus or other surface of higher genus) that cannot 
be lifted by any local perturbation. This degeneracy is topological in nature and it 
is intimately related to quantum number fractionalization. The robustness of the 
topological degeneracy against local noise due to the environment is at the core of 
topological quantum computation Q. 

Strong correlations that lead to these exotic quantum spectral properties can also 
impose kinetic constraints similar to those studied in the context of glass formers. 
Here we presents concrete examples of systems that have local Hamiltonians with 
no quenched disorder, exactly solvable spectra and topologically ordered quantum 
ground states. In these examples, not only do the classical thermal barriers grow 
with decreasing temperature, but also the widths of the quantum tunneling barriers 
do, in such a way that quantum processes are even more severely suppressed than 
classical ones at low temperatures. The origin of this behavior is the fact that any 
bath couples locally to the physical degrees of freedom of the system and it can 
only flip large objects through virtual processes of large order in the system-bath 
coupling. For these quantum systems, no physical bath can provide a mechanism 
for relaxation that is not dynamically arrested at low temperatures. 

We first discuss a two-dimensional (2D) quantum system with strong glass-like 
relaxation times when in contact with a restricted class of thermal baths; this 
example is used to clarify the issue of how a given Hamiltonian requires a minimum 
number of degrees of freedom that the bath must locally control for the system to 
be able to equilibrate. The second example is a three-dimensional quantum system 
with strong glass-like relaxation times for any class of baths that couple locally to 
physical degrees of freedom of the system. The third example is a three-dimensional 
quantum system with fragile glass-like relaxation times. 
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Figure 1. Square lattice, with spin operators defined on the sites. The square lattice is bipartite, and the 
two sets of points A,B are shown in red solid dots and blue open circles. A diamond contains 4 vertices in an 
elementary plaquette, and the diamonds can also be divided into two sets (forming a red/blue checkerboard) 
according to which sublattice their topmost vertices belong to. Four-spin operators are defined on each 
plaquette using the cr" and components of the spin, as described in the text. The green dots correspond 
to "defects" that are generated by applying a ct^ to the site encircled. 



3.1. Warmup: 2D example 

The first model is constructed on a two-dimensional (2D) square lattice, shown in 
Fig. [TJ Each site can be labeled by i,j G Z that index a site in the Bravais lattice 
spanned by the primitive vectors ai = (1, l)/\/2 and a2 = (—1, l)/\/2. To shorten 
the notation, define a superindex I = At every lattice site / one defines 

quantum spin S = 1/2 operators af, fij, and The square lattice is bipartite: 
it contains two sets of sites, which we label A and B, and which are shown in red 
and blue color in Fig. [TJ 

Let us define the quantum Hamiltonian in terms of the spins crj. Here we follow 
an approach similar to that of Kitaev, who constructed in a beautiful paper model 
quantum Hamiltonians that are exactly solvable In those models, the spins 
resided on links in planar lattices, but it is possible to carry out similar construc- 
tions with spins defined on vertices [iH], as it is done here. Later on we show how 
the construction with spins on vertices can be generalized to 3D lattices. 

Define a diamond cell Pj as the set of four lattice sites in an elementary plaquette 
with site / at its top. The four vertices are indexed by Jn{I)i for n = 1, . . . , 4, with 
one of the vertices Ji (/) = /. The four labels are assigned in such as way that the 
pairs { Ji, J3}, { J2, «^4} are diagonally opposite sites from one another. Explicitly, 
Ji(/) = / = J2{I) = (i - 1, j), = (i - l,i - 1), and J4(/) = {i,3 - 1). 

It is simple to see that the total number of diamonds equals the number of spins: 
each lattice site / is the top vertex of a single diamond. The one-to-one relation 
between a site / and the diamond Pi allows us to partition diamonds into two sets 
A and B (and color the corresponding diamonds red and blue, as shown in Fig. [1]). 

Now define the operators O/ as 

= ^J.(/) ^J3(/) ^J4(/) • (1) 

These operators commute, [Oj^Op] = for all pairs /, It is simple to see how: 
two diamonds Pi and Pp can share 0,1, or, 2 spins (/ 7^ /'). If they share spins, 
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they trivially commute. If they share 1 spin, the component (x,y or z) of cr for that 
shared spin coincides for both Oj and Oj' (the two diamonds touch along one of 
their diagonals). If they share 2 spins, the components of cr used in the definition 
of Oj and Oj' are different for both spins, there is a minus sign from commuting 
the X and y components of the spin operators from each of the shared spins, and 
the two minus signs cancel each other. 
Consider the system Hamiltonian 

H = -h^Oi, (2) 

which is trivially written in terms of the Oj operators, but complicated in terms of 
the original spins crj. Because the Oj commute, the eigenvalues of the Hamiltonian 
can be labeled by the list of eigenvalues {Oj} of all the Oj. Note that Oj = I, and 
each O/ = ±1. In particular, the ground state corresponds to O/ = 1 for all /. 

Because the number of spins equals the number N of sites, one may naively expect 
that the list {Oj = ±1} exhausts the 2^ states in the Hilbert space, spanned by 
{cr| = ±1}. However, there are constraints that the Oj satisfy when the system is 
subject to periodic boundary conditions (compactified to a torus). One can show 
that 

llOj=llOj = I. (3) 

leA leB 

There are two constraints; therefore there are only 2^~^ independent {Oj = ±1}. 
This implies, in particular, that there is a ground state degeneracy of 2^ = 4. 
(Notice that the ground state degeneracy is not associated with a symmetry. In 
particular, it is easy to show that {cr^'^'^) = 0.) This is a topological degeneracy 
and the eigenvalues of a set of two non-local (winding or topological) operators 7i,2 
are needed to distinguish between the 4 degenerate ground states. 

The operators 71,2 can be constructed as follows. Let Vi = {I \ i + j = 1} he a 
set of points along a horizontal line. Notice that sites on a line belong either all to 
sublattice A or all to sublattice B, for example Vi C A and V2 C B. Define 

Ti = n ^^=11 (4) 

i€Vi ier2 

It is simple to check that [7i,2, Oj] = for all / and the two operators 7i,2 trivially 
commute. Hence the two eigenvalues Ti^2 = il of 71,2 can distinguish the 4 de- 
generate ground states. This degeneracy has a topological origin and it scales with 
the genus of the surface. 

The spectrum of the Hamiltonian Eq. ^ is that of a trivial set of — 2 free 
spins, determined by the list of eigenvalues {0/ = ±1} of all the Oj, subject to 
the condition Eq. ([3]): E^Q^y = —h^^jOj. Excitations above the ground state 
(Oj = 1 for all /) are "defects" where O/ = — 1 in certain sites /. Because of the 
constraints in Eq. ([3j), the defects appear only in pairs. These defects have non- 
trivial quantum statistics: they are Abelian anyons with statistical angle different 
from that of fermions or bosons . 

The equilibrium partition function (within a topological sector) is given by 
Z = X]{Oj=±i} ^^^^^ thermal equilibrium at temperature T, the thermal 

average (Oj) = tanh^, and the concentration or density of 0/ = —1 defects is 
c = ^ (1 — tanh Y^)- Notice that we have encountered a situation analogous to clas- 
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sical spin facilitated models [i*] , in particular the plaquette models displaying glassy 
dynamics [H, iB-S]: the thermodynamics is trivial in terms of non- interacting de- 
fect variables. 

What about the dynamics of our quantum model? Although the spectrum of the 
model we are discussing is the same as that of free spins in a uniform magnetic 
field h, the variables Oj for different diamonds / cannot be independently changed. 
In fact, the operators Oj involve four spins, shared by neighboring operators Oj'. 
Under the action of any local spin operator, one cannot change the eigenvalue of 
Oj without changing the eigenvalues Oj' of its neighbors. What are the allowed, 
physical dynamical evolution rules for this quantum system? Can these dynamical 
rules lead to equilibration? 

In order to endow the system with some physical dynamics, we couple the orig- 
inal physical spins to individual baths at temperature T. Here we do not consider 
"Turkish" baths of multiple spins; nonetheless, as long as the groups of spins shar- 
ing a bath are locally delimited in space, the results obtained below should remain 
qualitatively unchanged. Moreover, allowing the bath to communicate informa- 
tion through long-ranged couplings (via phonons, for instance) will not change the 
results, as long as it operates on delimited regions of space. 

When the original individual spins are coupled to their baths, "flips" of the states 
of multiple Oj sharing a given spin take place. Therefore, our model has a trivial 
spectrum but a highly correlated dynamics. It is this correlated dynamics that 
gives rise to non-trivial non-equilibrium behaviour. 

We introduce the bath deg rees of freedom as in the Feynman- Vernon influ- 
ence functional approach [19l | or Caldeira-Leggett dissipative quantum mechan- 



ics formulation [201 . l2l[ |. by letting the Hamiltonian of the system plus bath be 



% = H + i/bath + ^^spin+bath, where H is defined in Eq. ([2]), and 

/■oo 

i^bath = dx [n?(i, x)]2 + [d.^fit, x)? (5a) 

4pin/bath = Y.9aCTf nf{t, 0) . (5b) 
I,a 



The three components (a = 1,2,3) of the conjugate vector fields and 11/ obey 
the equal-time commutation relation [<I>"(i, x), IIj {t,x')] = i 5ij 6aa' ^ix — x'). 

Notice that, for each site /, the bath-spin system can be viewed as an extended 
bosonic string that couples to a spin at the boundary x = 0. The coupling ampli- 
tudes are ga- One can in general choose anisotropic couplings, but the most general 
bath should contain all of 51,2,3- In the quantum model, acting on a site /' G Pj 
with one of ap , crj, , or a}, flips or not the eigenvalue Oj depending on whether 
i^x,y,z _ cFp^'^, respectively. 

If integrated out, the bath degrees of freedom away from the boundary x = 
lead to a non-local in time action and to dissipation effects. Instead of working 
with the dissipative action, let us follow the time evolution of the system plus bath 
and look at the possible evolution pathways of the quantum mechanical amplitudes 
of the system plus bath degrees of freedom. After evolution by time t from some 
initial state, the system is in a quantum mechanical superposition 

l*)= E r{o.} l{0/})^|T{o,}), (6) 

{Oj=±l} 

where |T|o^|) is a state in the bath Hilbert space with norm one. (Here we focus 
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on states in a single topological sector, although mixing sectors can be done by 
including the eigenvalues of the topological operators T; mixing is exponentially 
suppressed as the system size increases) . The fact that the bath degrees of freedom 
couple to single quantum spins crj enters in the problem through the permitted 
channels for transferring amplitudes among the TjO/}- 

The processes that transfer amplitude among the ^{Oi} correspond to different 
orders in perturbation theory on the system-bath coupling g^. There is also a 
thermal probability factor coming from the bath and that depends on the difference 
between the initial and final energy E^Q^y = —h^^Oi of the system. One class 
of paths is a sequential passage over states connected through order Qa processes; 
this is a "semi-classical" type trajectory. 

Within this restricted class of processes, we can make a connection to the classi- 
cal plaquette models, particularly the 4-spin square plaquette model whose glassy 
properties have been studied in Refs. H, 12, 23]. This classical model is obtained 



by defining diamond variables r/ in place of the Oj as in Eq. ([T|) using only, say, 
the z-component cr^ for all four sites of the diamonds. Flipping an individual spin 
changes signs to all four r/ variables surrounding the spin. This multi-defect type 
dynamics makes it difficult for the system to relax to equilibrium. For example, 
if the temperature is lowered, in order to decrease the defect density, either four 
defects come together and annihilate (4 — )• decay), or three defects become one 
(3 — )• 1 decay). However, the defects are not free to diffuse and come together. 
In order to move an isolated defect, it must first decay into three defects (1 — )• 3 
production) , then a pair can diffuse freely, and recombine with another defect else- 
where through a 3 — )• 1 decay process. The final result is that the original defect 
(as well as the other defect elsewhere) moves by one lattice spacing and the total 
number of defects does not change (it first increases by 2 and then decreases by 2). 
Because of the initial 1 — )• 3 production process, there is an energy barrier of 2h to 
overcome. This activation barrier leads to recombination/equilibration times 

tseq. ~ exp(2/i/T) (7) 



that grow as temperature is lowered in an Arrhenius fashion [Ji 

In our model, the same situation is recovered if the coupling to the bath involves 
only cr^ components of the spins (i.e., 5i = 52 = 0). The Hamiltonian encompasses 
only (T^ and components and the action of a operator on a spin changes 
the eigenvalues of all four surrounding O/. One may then wonder how quantum 
tunnelling processes modify the behaviour of the system as compared to its classical 
counterpart. Recall that defect annihilation occurs only through virtual processes 
in which the number of defects is strictly larger in the intermediate (virtual) steps. 
At temperature T, the typical defect separation is ^ = c~^/^ ~ e^/^"^. Bringing 
them together requires tunnelling processes at least of order ^ in perturbation 
theory, which have an amplitude of order [g/h)^ (notice the energy denominator 
h). This amplitude leads to recombination/equilibration times 



ttun. ~ exp 



ln(V5)e'^/2T ^ 



which grow extremely fast as the temperature is lowered. What we learn from this 
simple estimation is that quantum tunneling is less effective than classical sequen- 
tial processes in thermalizing the system. This is counterintuitive to the notion 
that at low temperatures quantum tunneling under energy barriers remains an 
open process while classical mechanisms are suppressed due to thermal activation 
costs. The reason for this particular quantum freezing is simple: equilibration pro- 
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ceeds through activated defect recombination; as the density of defects decreases 
at low temperatures, the barrier widths increase and debihtate quantum tunnehng. 
In passing, we note that in a finite system of size L, one must replace ^ by L in the 
estimation of the recombination/equilibration times, ttun. ~ exp [ln(/i/(jr) L]; this 
time scale is also of the order of that for tunneling between two topological ground 
states in a finite system of size L 

Let us return to the issue of which component of spin enters in the coupling to the 
bath. The minimal bath coupling involves the z-component of the spin operator, 
for this component does not commute with any of the Oj that contain that given 
spin. However, there is no a priori reason why the coupling should be restricted 
to z-components only. If the bath couples to x- and y-components as well, defects 
can diffuse freely without barriers (via 1 — )• 1 processes), eliminating the need for 
1 — )• 3 defect production processes in order to bring them together and annihilate 
them. Simple defect diffusion leads to a diffusive equilibration time t^g ~ for a 
system of size L (polynomial in L with constant exponent). 

Such dependence on the details of the bath coupling is removed in the 3D models 
we discuss next, one of which has even slower equilibration, as in fragile glasses. 



3.2. 3D quantum strong glass 

In the previous example, particle diffusion finds its origin in the fact that the x- or 
z-components of a spin commute with 2 out of the 4 surrounding diamonds. There- 
fore, defects can be created in pairs, not quadruplets, and single defect diffusion 
can take place through annihilation of a defect and creation of a neighbouring one 
via pair flip. In the example that follows (see also Refs. [1, US]), six defect cells 
share each single spin, in such a way that acting with any component of the spin 
operators flips 4 cells and defect diffusion cannot occur. 

The model displaying strong like glassiness is constructed on a three-dimensional 
(3D) face-centered cubic (fee) Bravais lattice, spanned by the primitive vectors 
ai = (l,l,0)/\/2, a2 = (0,l,l)/\/2, and ag = (1,0, l)/\/2. Each site can be 
indexed by i,j, A; G Z, and to shorten the notation, define a superindex / = {i,j, k). 
At every lattice site / one defines quantum spin 5 = 1/2 operators af, a^, and o"|. 

The fee lattice hosts sets of octahedra: the simplest one to visualize is the one 
assembled from the centers of the six faces of a cubic cell, as shown in Fig. [2l In 
addition to this simple set, there are three more sets of octahedra that can be 
assembled from sites both on faces and on corners of the cubic cells, totalling 4 
such sets, which we label A, B, C and D. 

It is simple to see that the total number of octahedra equals the number of spins: 
each lattice site I is the topmost vertex of a single octahedron. Define then Pj as 
the set of six lattice points in the octahedron with site / at its top. The six vertices 
are indexed by Jn(/), for n = 1, . . . , 6, with one of the vertices Ji{I) = I. The six 
labels are assigned in such as way that the pairs {Ji, J4}, {J2,Jr>}-, {<^3)>^6} are 
diagonally opposite sites from one another. This number labeling is illustrated for 
a single octahedron in Fig. [2j From the one-to-one relation between a site I and 
the octahedra Pi, we can also partition the lattice sites into the four sets A, B, C 
and D of octahedra. 

Now define the operators O/ as 

= ^Jl(/) ^J.(7) ^J3(/) '^^4(/) ^J5(/) ^J6(/) • 

These operators commute [Oi,Oi'] = for all pairs Indeed, two distinct 

octahedra Pi and Pp can either share 0,1, or, 2 spins. If they share spins, they 
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Figure 2. Cubic cell of an fee lattice. The centers of the six faces form an octahedron, with its sites 
labeled from 1 (topmost) to 6. In addition to the set of octahedra formed by the face centered sites, there 
are three more sets of octahedra that can be assembled from sites both on faces and on corners of the cubic 
cells, totalling 4 such sets. Six-spin operators are defined on these octahedra using the a^'^''^ components 
of spin on each vertex as described in the text. 



trivially commute. If they share 1 spin, the component (x,y or z) of a for that 
shared spin coincides for both Oi and Oi' (the two octahedra touch along one of 
their diagonals). If they share 2 spins, the components cr used in the definition of 
Oi and O/' are different for both spins. Thus, there is a minus sign from commuting 
the spin operators from each of the shared spins, and the two minus signs cancel. 

Consider the system Hamiltonian as in Eq. ([2|) , which is trivially written in terms 
of the Oi operators, but complicated in terms of the original spins <t/. Because 
the Oi commute, the eigenvalues of the Hamiltonian can be labeled by the list 
of eigenvalues {O/} of all the O/. Notice that Oj=I and so each Oj = ±1. In 
particular, the ground state corresponds to O/ = 1 for all /. 

Because the number of spins equals the number N of sites, one may naively expect 
that the list {Oj = ±1} exhausts the 2^ states in the Hilbert space. However, there 
are constraints that the Oj satisfy when the system is subject to periodic boundary 
conditions (compactified) . One can show that 

l[Oi=llOj = llOj=llOj = l. (10) 

leA I(^B I&C I&D 

There are four constraints; therefore there are only 2^~^ independent {O/ = ±1} 
and there is a ground state degeneracy of 2^ = 16 (for a thorough discussion of 
this degeneracy and its dependence on boundary conditions, see Ref. fill). This is 
a topological degeneracy and the number of ground states scales with the genus of 
the manifold the system is defined on. The eigenvalues of a set of four non-local 
(topological) operators 71,2,3,4 are needed to distinguish between the 16 degenerate 
ground states. 

The operators 71,2,3,4 can be constructed as follows. Let Vi = {I \ j + k = 1} he a. 
set of points along a horizontal plane. Notice that each plane contains sites in only 
two of the four sublattices A,B,C,D. For example Vi C AU B and V2 C CU D. 
Define 

Ti= H T2= H %= H %= H (11) 

/e-PiHA leViiiB ieV2nc i&V-inD 

It is simple to check that [71,2,3,4, C/] = for all /, and the 71,2,3,4 trivially com- 
mute among themselves. Hence the four eigenvalues Ti, 2,3,4 = ±1 of 71,2,3,4 can 
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Figure 3. Sites of an hep (hexagonal close-packed) lattice, (a) The hep lattice is comprised of two 
interpenetrating hexagonal lattices, which can be alternatively seen as stacked triangular lattices, shown 
in red and blue. Prisms are defined as sets of five sites, two of which belong to one sublattice (top and 
bottom of the prism), and three of which belong to the other and form a triangle that lies in the layer 
in between the top and bottom sites of the prism. Five-spin interactions are defined on each prism as 
explained in the text, (b) Top view of the hep lattice, which shows that the blue sites stack on top of 
the red upward pointing triangles, and the red sites stack on top of the downward pointing blue triangles, 
(c) Two prisms with topmost sites belonging to different sublattices can share a common edge and the 
five-spin operators defined on the two prisms commute because minus signs from commuting the and 
cr^ components appear twice, once for each shared site, and cancel. 

distinguish the 16 degenerate ground states. 

In this model one can verify that, whichever component of spin enters in the 
couphng to the bath, it is impossible to flip only a pair of defects and thus there 
is no mechanism for defect diffusion. The reason is that any site is shared by 6 
octahedra, and the operators Oj for these cells are such that one can divide the 
6 into 3 groups of 2 octahedra that will have in their definitions, respectively, the 
X, y, and z component of spin operators at the shared site. Acting with either of 
the three components of the spin operator on this shared site will flip at least four 
defects. 

Quantum glassiness, i.e., the behaviour in Eq. ([8|) leading to Tq ~ exp(L), is thus 
protected against any local bath. On the other hand, as soon as we leave the T = 
limit, the system size dependence of the relaxation time scales changes radically. 
Indeed, the energy barrier to defect diffusion is only a 'one-step' process. As soon as 
1 — 3 decays are allowed to take place thermally, any two of the three new defects 
can diffuse freely across the system. The time scales for diffusion are simply rescaled 
by the factor tseq. ~ exp(2/i/T), Eq. ([7]), but the system size dependence remain 
only Tc ~ poly(L). 

3.3. 3D quantum fragile glass 

The model displaying fragile-like glassiness is constructed on a three-dimensional 
(3D) hexagonal close-packed lattice, shown in Fig. O The lattice can be viewed as 
two interpenetrating simple hexagonal Bravais lattices displaced from one another 
by ai/3 + aa/S + 03/2, where ai = (1,0,0), aa = (l,\/3,0)/2, and ag = (0,0,1) 
are the primitive vectors of the simple hexagonal lattice. The sites belonging to the 
two intercalating lattices are shown in red and blue color in Fig. [31 Each site can 
be labeled by k £ Z that index a site in the Bravais lattice spanned by ai,2,3! 
plus q = 0,1 that indexes each of the two sublattices ~ to shorten the notation. 
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define a superindex / = k; q). At every lattice site / one defines quantum spin 
S = 1/2 operators o"^, fij, and ct|. 

Define now a prism cell Pi that contains five lattice sites Jn{I), for n = 1, . . . , 5 
as follows. For a given lattice site /, the prism Pj contains the site Ji{I) = I, 
which belongs to one sublattice of the hexagonal close-packed lattice, the three 
sites that belong to the other sublattice and that form an elementary triangle 
(sites J2, J3, J4) just below the site /, and finally the site J^il) just below that 
triangle, which belongs to the same sublattice as site /. [In terms of the lattice 
indices, J5(/) = (i, j, k — l;q) if / = (i, j, k; q).] An example of two prisms is shown 
in Fig. [3]3. Notice that the two prisms shown share a common edge, and that their 
tops belong to distinct (red and blue) sublattices. The total number of prisms 
equals the number of spins: each lattice site / is the top vertex of a single prism. 

Now define the operators Oj as 



= ^Ji(/) "^Mi) ""MI) ""MI) ^J5(/) • (12) 



The operators commute, [Oj,Oi'] = 0, for all pairs /, Indeed, if /, /' belong 
to the same sublattice and the prisms Pj,Pj' share a vertex, then they trivially 
commute as they both involve the same component (x or z) of the spin operators 
a at the shared site. If they belong to distinct sublattices, they either share spins 
or an edge with 2 spins, as shown in Fig. [3j If they share 2 spins, the minus signs 
from commuting the x and z components of spin in each of the shared sites appear 
an even number of times. 

Consider the system Hamiltonian as in Eq. ([2]) , which is trivially written in terms 
of the C/ operators, but complicated in terms of the original spins aj. Because 
the Oj commute, the eigenvalues of the Hamiltonian can be labeled by the list 
of eigenvalues {0/} of all the Oj. Notice that Oj = I, and so each 0/ = ±1. In 
particular, the ground state corresponds to 0/ = 1 for all /. 

Because the number of spins equals the number N of sites, one may naively expect 
that the list {Oj = ±1} exhausts the 2^ states in the Hilbert space. However, there 
are constraints that the Oj satisfy when the system is subject to periodic boundary 
conditions (compactified). Each of the two sublattices (g = 0, 1) of the hexagonal 
close-packed structure can be further subdivided into Ag, Bg or Cg according to the 
three sublattices of the tripartite triangular stacks of the simple hexagonal lattice. 
(All in all, there are six sublattices Aq^i,Bq^i and Co,i.) One can show that 



n oj= n n ^^=^- (i^) 

leAgUBg I&BgUCg /GC,UA, 



These are six constraints, but only four are independent, because the product of the 
three products in Eq. (I13p for the same q is trivially the identity. Therefore there 
are only 2^~^ independent {Oj = ±1}. This implies, in particular, that there is a 
ground state degeneracy of 2^ = 16 which is topological in nature and the number 
of ground states scales with the genus of the manifold the system is defined on. 
The eigenvalues of a set of four non-local (topological) operators 71,2,3,4 are needed 
to distinguish between the 16 degenerate ground states. 
The operators 71,2,3,4 can be constructed as follows. Let the plane Vk,q be the 
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set containing sites with fixed k and q. Let 

ri= n r2= n (14^) 

ieVi,on(AoUBo) /ePi,on(BoUCo) 

r3= n ^1 '7'4= n ^i- (14b) 

/ePi.in{AiUBi) /ePi,in(BiUCi) 

It is simple to check that [71,2,3,4,0/] = for all I, and the 71,2,3,4 trivially 
commute among themselves. Hence the four eigenvalues Ti,2,3,4 = ±1 of 71,2,3,4 can 
distinguish the 16 degenerate ground states. 
There are relations between this 3D model and a 2D classical triangular plaquette 



model which has glassy behavior [IJ, ll6|, UJ^ ■ The 2D triangular plaquette model 



has Ising spin variables defined on the sites of a triangular lattice and a 3-spin 
interaction which is the product of the Ising variables on the downward pointing 
triangular plaquettes only. Plaquette Ising variables (the 3-spin products) are de- 
fined at the centre of the downward triangles, which behave thermodynamically as 
free Ising spins in a magnetic field. However, the dynamics is rather non-trivial in 
terms of the plaquette variables, for flipping an original spin corresponds to flipping 
all three surrounding plaquettes. 

In our 3D model, each quantum spin crj is shared by 5 prisms: 3 whose centers are 
on the same plane, and 2 whose centers are immediately above and below site /. If 
the system-nath coupling contains the spin component, all 5 prisms are flipped. 
The and components flip either the eigenvalues of the 3 prisms on the plane 
or the 2 prisms on the vertical direction, respectively. Flipping the eigenvalues of 2 
prisms in the vertical direction leads to defect diffusion, but only in that direction. 

To connect our 3D quantum model to the triangular plaquette model, consider a 
compactified slab (periodic boundary conditions) in the third dimension (parallel 
to as), with M layers. Because of the periodic boundary conditions, the odd-even 
parity of the defect number is conserved along vertically stacked prisms regardless 
of the system-bath spin-flip operator, a^, a^, or a^. The defect number parity can 
be captured by defining the following operator (recall I = k; q)): 



It is also useful to define a similar product over the third dimension for the original 
spins: 



These "slab" operators allows us to concentrate on subspaces of the Hilbert space 
with a given set of Ti^-^q instead of the states with given Oij^k;q- There are dynamical 
processes that transfer quantum mechanical amplitudes within and between these 
subspaces labeled by Tij-g] we can argue that the system is glassy by simply looking 
at the processes that transfer amplitude between the subspaces. 

The variables Tij-g and Sij-g can effectively be used to relate our quantum model 
to two 2D systems (g = or red, and g = 1 or blue) defined on sites labeled by 
q) of two distinct triangular lattices. The variables Sj j-g can be related to the 
original spin variables in the models of Refs. [13, [H, Ea]. In particular, one can 
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relate the and the Tij-^q using Newton's binomial coefficients through 



n,m\q\ 



(17) 



Using the Pascal triangle relation 
indeed the defect variables correspond to 



+ 



one can show that 



(18) 



We can focus again on two classes of processes: sequential passage over states con- 
nected through order ga processes ("semi-classical" type trajectories), and quantum 
tunneling process. The a naly sis of the sequential processes is similar to that of the 
classical 2D models lid . ITtI and goes as follows. Single flips of s correspond 
to concomitant flips of three r defects. Defects can only be annihilated in triplets. 
The defects are not free to diffuse and come together; instead, they move through 
the production of more defects. For example, a defect can decay into two more 
defects, by flipping one s variable. Now, in order to bring three defects separated 
by a distance ^ together, one has to go through intermediate steps with a large 
number of defects that are created. There is a hierarchical organization of these 
intermediate processes; equilateral triangles of size ^ = 2^ require the creation of 
i extra intermediate defects. Hence there is an energy barrier of order ih to be 
overcome. For a typical equilibrium separation ^ = c~^/^ ~ e'*/^^, the barriers to 
be overcome in the equilibration processes are of order h'^/{T2ln2). Hence, the 
equilibration time scales as 



''seq. 



exp 



21n2 



(19) 



a much slower relaxation than the Arrhenius one, Eq. ([7]), for the model discussed 
in Section 13. 2[ 

Through quantum tunneling, defect annihilation can occur via virtual processes 
in which the number of defects is strictly larger in the intermediate (virtual) steps. 
The order in perturbation theory in g grows very fast with defect separation. An 
example is shown in Fig. HJ basically, to annihilate three defects at the corners of an 
equilateral triangle of size ^ = 2^, one must flip 3^ original spins laying on a mold 
defined by a Sierpinski gasket. (Notice that here the hierarchy is built staring from 
the microscopic scale.) In perturbation theory the quantum recombination process 
has an amplitude of order (g/h)^ , which leads to recombination/equilibration times 



''tun. 



exp 



ln(^ 



exp 



ln3 /i 
21n2r 



(20) 



which grows extremely fast as the temperature is lowered. Again, we learn from this 
simple estimation that quantum tunneling is less effective than classical sequential 
processes in thermalizing the system. 
In a finite size system, the limiting time scale will be set by ^ ~ L, whereby 



exp 



Ll°&(3) In ( I 
h 



(21) 



leading to a quantum mechanical exponential dependence on system size in contrast 
with the polynomial scaling of the classical (activated) time scales that become 
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Figure 4. To annihilate three defects (shown in green) at the corners of an equilateral triangle, one must 
flip the spins in a "fractal" membrane (containing sites shown in red) that stretches between the defects. 
For a triangle of size 2*, there are 3* sites in the membrane. The annihilation of the three defects through 
quantum tunneling is a virtual process of order the number of sites that are involved (red sites). Hence, the 
amplitude for the quantum tunneling process vanishes exponentially with the "volume" of the membrane. 



available at any finite temperature, 



Tc ~ exp 



/ilog2(L) 



In 2 



(22) 



as in the long time scale limit ^ — )• L, the classical energy barriers ih tend to the 
value hlog2{L). 



Note on perturbative approaches 



It is interesting to comment briefly on more generic approaches to study equili- 
bration in quantum systems coupled to a bath, specifically in regards to quantum 
glasses with exponential time scales. 

Quantum systems in contact with a reservoir are characterized by mixed en- 
sembles of states that are best described using von Neumann's density matrix 
formalism 25|, [26]. At equilibrium, the density operator is given by p = e~^^ /Z in 
the canonical ensemble, where H is the Hamiltonian of the system, /3 = 1/T is the 
inverse temperature of the reservoir (and of the system of interest, at equilibrium), 
and Z = tr(e~^^) is the canonical partition function. The expectation value of any 
operator O in this mixed ensemble is given by tr(/50). Therefore, the characterisa- 
tion of the properties of a quantum system at equilibrium is essentially a spectral 
problem. Describing the low temperature properties of a system, for example, re- 
quires the understanding of the ground state, its symmetries (or lack thereof) and 
its quantum orders [27|, and its low lying excitations. In particular, at zero tem- 
perature, the system goes to its ground state. In the case of a quantum many-body 
system, if parameters in the Hamiltonian H are tuned, the ground state can change 
its symmetries or its quantum orders through quantum phase transitions. 

The time evolution of a system s in contact with a reservoir b and total Hamil- 
tonian 



H = Hs + Hi, + V = Ho + V 



(23) 
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is described by the Liouville-von Neumann equation 

dp{t) 



dt 



i[H,p{t)], (24) 



where p{t) = exp{iHt) p exp{—iHt) is the time-dependent density matrix of the 
system plus bath, and we use the convention h = 1. (Without loss of generality, 
we assume that H does not depend directly on time.) It is customary to recast 
Eq. (j24p in the interaction picture format, 

pi{t) = e-'^"^ e*^* p e-*^* e*^"* (25) 
Vi{t) = e-'^°^ V e'""^ (26) 



whereby 

dpi{t) 
dt 



-i[Vi{t),pi{t)]. (27) 



Note that /0/(0) = /3(0) = p. If we are interested in the evolution of the system 
alone, we need to trace over the degrees of freedom of the reservoir, 

pf{t)^tn[pm (28) 

^^2^ = -itr, [!/,(*), (29) 

The right hand side of Eq. (I29p results in a linear super-operator acting on p^j^ (t) . 
By comparison with its classical counterpart, i.e., the transition matrix in a Master 
equation, the eigenvalues of the super-operator determine the equilibration rates of 
the quantum system coupled to the bath. Following Ref. [11], one could then sym- 
metrise the super-operator and reinterpret it as a fictitious quantum 'Hamiltonian' 
acting on the space of density matrices p. The appearance of diverging time scales 
in the quantum system would then correspond to a quantum phase transition in 
the associated 'Hamiltonian'. Note that at the quantum level, the correspondence 
between the original Hamiltonian and the fictitious one occurs interestingly in the 
same number of dimensions for both systems. 

Unfortunately, taking the trace of the right hand side of Eq. (j29l) to find the super- 



operator acting on p\^\t) is in general a tall order. A common approach consists of 
expanding the commutator [Vi(t), pj{t)] perturbatively in the interaction between 
system and reservoir, leading to the recursive equation 

+ E(-^)"+' fdn / V2 . . . r^drn [Viit), [Viin), . . . [ViiTn),pim] . . .] 
Jo Jo Jo 

= -i[Vj{t),pjm- fdT[Vj{t),[Vj{T),pjm] + (30) 

Jo 

Carrying out the trace becomes now feasible under the conventional assumption 
that the initial density matrix factorises into system and bath, pi{0) = p^'^^ <^ p^^\ 
However, one needs to consider that each order in the expansion enables tunnelling 
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processes in the system-bath interaction V up to that very same order. As discussed 
in Sec. [2] and as seen in the concrete examples in Sec. [3l the relaxation processes 
typical of quantum glasses require tunnelling under barriers whose width scales with 
system size. Such processes would be all together forbidden at any finite order in the 
perturbative expansion. One could view this as a signature of quantum glassiness: 
the appearance of disconnected sectors in the quantum dynamics of the system at 
any finite order in perturbation theory. 

5. Conclusions 

There are plenty of systems in nature that recalcitrantly avoid equilibration with 
the environment. When their time scales become impractically large to measure, 
these systems are conventionally referred to as glasses. 

In order to investigate how this phenomenon comes about, it is reasonable to 
distinguish between different classes of slow dynamics according to the way that 
their equilibration times scale with their physical size. One can draw a parallel 
with computational complexity. When the system does not manage to equilibrate 
in poly(L) time, it is because nature's algorithm for (local) dynamics is not efficient. 
Systems with exp(L) time scales are computationally hard given nature's resources. 
[We note in passing that, just as in the case of the definition of computational 
complexity classes, one must allow for considerations of practical importance: time 
scales can diverge as a function of temperature while remaining polynomial in L.] 

What happens when temperature is lowered all the way down to T = 0? Because 
absolute zero temperature freezes out all classical activation processes, quantum 
tunneling is all that is left to dynamics. The relevant issue then becomes which 
kind of tunneling barriers remain to surpass so as to reach the ground state in 
the presence of a zero temperature bath. It is the nature of these barriers that 
separates systems with relaxation times that are exponential vs. polynomial in 
L, and therefore separates systems that are quantum glasses (hard, for nature's 
algorithm) from those that are not. 

In this paper we discussed these issues and constructed explicit examples of sys- 
tems with exponential time scales at zero temperature. These systems are devoid 
of disorder or local symmetry breaking, but rather exhibit a non-local form of or- 
der known as topological order. As such, they ought to be considered examples 
of topological quantum glasses. We showed how our local quantum Hamiltonians 
resist equilibration with a thermal bath and glassiness is "protected" against ther- 
mal equilibration with any bath, provided that it couples locally with the physical 
degrees of freedom of the system. 

It is often believed that quantum tunneling provides an escape route against 
classical dynamical slowdown caused by thermal energy barriers as the temperature 
is lowered. In the systems presented here, classical sequential processes are more 
effective than quantum tunneling processes in the limit T — )• 0. The reason for 
the freezing of quantum tunneling is that equilibration is through activated defect 
recombination; as the density of defects decreases at low temperatures, the barrier 
widths increase and eventually become as large as the size of the system. 

We also discussed briefly generic approaches to study equilibration in quantum 
systems in contact with a reservoir, starting from the Liouville-von Neumann equa- 
tion. Following this route one can in principle arrive at the quantum mechanical 
counterpart to the classical Master equation for the equilibration of probabili- 
ties. The role played by the transition matrix is here taken by a super-operator, 
whose spectrum controls the relaxation properties of the system. It is intriguing to 
speculate that a quantum-to-quantum mapping akin to the classical-to-quantum 
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correspondence discussed in Ref. jl^ would lead to the formulation of a fictitious 
quantum mechanical system in the same number of dimension of the original one, 
where quantum dynamical transitions would appear as static, zero temperature 
quantum phase transitions. However, constructing the super-operator is in general 
a tall order and at best one can do so order by order in perturbation theory. One 
should note however that truncating the expansion to any fixed order in perturba- 
tion theory on the system/bath coupling disconnects the space of states. In order 
to study quantum glassy systems, it is necessary to go to all the way to order L 
in the perturbation expansion, accessing matrix elements that are exponentially 
small, and therefore time scales that are exponentially large in L. 
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